Quantum melting of the quasi-two-dimensional vortex lattice in k— (ET)2Cu(NCS)2 
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We report torque magnetization measurements in regions of the mixed state phase diagram (B ~ 
/ioHc2 and Tc/10^) of the organic superconductor k— (ET)2Cu(NCS)2 where quantum fluctuations 
are expected to dominate thermal effects. Over most of the field range below the irreversibility line 
(Birr), magneto-thermal instabilities are observed in the form of flux jumps. The abrupt cessation 
of these instabilities just below Bin- indicates a quantum melting transition from a quasi-two- 
dimensional vortex lattice phase to a quantum liquid phase. 

PACS numbers: 74.70.Kn, 74.60. Ec, 74.60. Ge, 74.25.Ha 



Although the properties of vortices in layered type— II 
superconductors have been studied vigorously for over 
15 years, this subject continues to provide a rich and 
varied field for investigation. The mixed superconduct- 
ing state represents a wonderful playground for study- 
ing general phase transformations associated with vortex 
matter [Q. In particular, melting of the Abrikosov vortex 
lattice into a liquid phase has drawn considerable inter- 
est, especially in extreme type— II superconductors such 
as the high temperature superconductors Q (HTS) and 
organic superconductors Moreover, identification of 
a melting transition driven by quantum (as opposed to 
thermal) fluctuations has drawn much attention, both 
theoretically as well as experimentally A thor- 
ough understanding of vortex physics is also essential for 
determining possible limiting behavior for technological 
applications utilizing superconductors. 

A model system for investigation of the structure 
and dynamics of vortices in layered, type— II, supercon- 
ductors is K-(ET)2Cu(NCS)2, where ET denotes bis- 
ethylenedithio-tetrathiafulvalene |6|. Like the HTS, this 
organic superconductor possesses a highly anisotropic 
layered structure with the superconducting ET planes 
separated by insulating anion layers The anisotropy 
parameter, 7, defined as the ratio of the interlayer (cur- 
rents II a- axis) and in-plane (currents || 6c-plane) pene- 
tration depths (7 = Aa/Ahc), is thought to be in the range 
50-200 [|, i.e. similar to that of BiaSrsCaCuzOg+d %■ 
Consequently, in comparison to conventional supercon- 
ductors, fluctuation effects may be expected to play an 
important role at low-temperatures. In contrast to the 
HTS, the organic superconductors are extremely clean, 
with very few crystal defects. Furthermore, because of 
the reduced T^ and Hc2 (T^ = 9 - 10 K ||,|,|l3 and 
AioHc2 « 5 T (30 - 35 T) for the field perpendicu- 

lar (parallel) to the superconducting layers), one is able 
to probe much more of the temperature/field parameter 
space within the superconducting state than is currently 



possible in the HTS. This, in turn, raises the interest- 
ing prospect of studying two dimensional vortex lattice 
melting over an extended temperature range (from Tc 
to Tc/10'^) — possibly into the (low T, high B) quantum 
regime. Magnetic relaxation measurements performed by 
Mota et al. |l^, have shown that the crossover temper- 
ature from thermal to quantum dominated fluctuations 
occurs at around 0.5 K for this material. 

It is well known from muon spin rotation (/iSR) mea- 
surements that a three-dimensional (3D) flux line lattice 
exists only at very low flelds (< 7niT<C /ioHca) in the 
title compound |p^; due to the large 7 value, the individ- 
ual quasi-two-dimensional (Q2D) vortex lattices in adja- 
cent layers become effectively decoupled above a roughly 
temperature independent dimensional crossover field Bcr 
7mT, see Fig. 1 and Ref. flj]). Our recent tempera- 
ture dependent investigations of the interlayer Josephson 
plasma resonance (JPR) in k— (ET)2Cu(NCS)2 provide 
conclusive proof that, in spite of the decoupling, long 
range Q2D order among vortices within individual layers 
persists over most of the region in the phase diagram be- 
tween Bcr and the irreversibility line (see shaded region 
in Fig. 1 and Refs ]l3| , |l^ ). Immediately below the ir- 
reversibility line, local magnetization measurements [^, 
as well as the JPR studies ||l^, indicate a transition in 
this Q2D vortex structure (see Fig. 1) — Q2D melting 
has been put forward as one possible explanation for this 
transition 13|. These findings serve as a motivation for 
investigating a wider range of the B,T phase diagram. 

High quality single crystals were grown using standard 
techniques Q . A single sample (approximate dimensions 
1.0 X 1.0 X 0.3 mm'^) was mounted on a capacitive can- 
tilever beam torque magnetometer which, in turn, was 
attached to a single axis rotator; 6 = 0° corresponds 
to the field parallel to the least conducting a-axis, while 
9 = 90° corresponds to the field parallel to the highly 
conducting &c-planes. The sample, cantilever and rota- 
tor were then loaded directly into the mixing chamber 
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of a top-loading ^He/^He dilution refrigerator situated 
within the bore of a 20 T superconducting magnet at the 
National High Magnetic Field Laboratory (NHMFL). For 
all measurements reported here, the magnetic field was 
swept at a constant rate of 0.5 T/min. Temperature de- 
pendent torque measurements were performed at 6 — 47° 
and 74°. Subsequent analysis of the angle dependence 
(from 6* = 0° to 90°) enabled us to scale the temperature 
dependence back to 6 = 0°, where the torque is zero in 
this setup. 

In Fig. 2, we plot the magnetization M, derived di- 
rectly from the torque r (M — r/Bsin 9), as a function of 
the applied magnetic field strength B, for angles between 
6 ^ 30° and 79° at approximately 10° intervals; the tem- 
perature is 25 mK. The overall shape of these curves is 
consistent with previous measurements The ob- 

served magnetization and associated hysteresis (arrows 
in Fig. 2 indicate the field sweep direction) are a con- 
sequence of the viscous flow of magnetic flux into (out 
of) the sample upon increasing (decreasing) the applied 
magnetic field strength. The hysteresis is greatest at low 
fields, and can be seen to disappear completely above a 
characteristic field B^^r (< MoHc2), above which the sam- 
ple behaves reversibly. Perhaps the most pronounced fea- 
tures of the data in Fig. 2 are the abrupt magnetization 
jumps. These "flux jumps" have been observed previ- 
ously in this and other materials ||l^,|l^,^, and are due 
to an avalanche behavior associated with the reorganiza- 
tion of magnetic flux as it enters the sample; a systematic 
analysis of this phenomenon can be found in Ref. jl^ . 

In the region between Bcr and Hirr, crystal defects 
collectively pin the Q2D vortex lattices in each layer. 
Thus, in a field swept experiment, there is a build up 
of flux near the sample surface. This creates a field 
gradient at the sample edge together with an associ- 
ated surface current given simply by Maxwell's equation: 
V X B = /ipJe, where Jc is the in-plane critical current 
density [g0|j2l|]. At extremely low temperatures, a ther- 
mal boundary (Kapitza) resistance isolates the sam- 
ple from the surrounding cryogen bath. Viscous trans- 
port of vortices across the sample edge (where the sur- 
face screening currents flow) causes local heating which, 
in turn, reduces the critical current density, leading to 
additional heating and so on.. In the absence of an 
effective thermal link to the surroundings, this can re- 
sult in runaway thermal instabilities which cause macro- 
scopic regions of the sample to become metallic |l8| , p3| . 
When this occurs, flux is able to flow rapidly into the 
crystal's interior, thereby negating any/most of the sur- 
face currents. The crystal then quickly cools and once 
again becomes superconducting with a slightly different 
metastable vortex arrangement; the process then starts 
anew. Indeed, Legrand et al. p^ , have shown that, for a 
field swept experiment, YBa2Cu307 samples experience 
large temperature spikes in conjunction with discontinu- 
ities in the magnetization. These jumps are then followed 



by a rapid relaxation back to the surrounding bath tem- 
perature, consistent with the model above. An analysis of 
the temperature dependence of the magnitudes of the ob- 
served magnetization jumps (ocT'^/^) in Fig. 2 supports 
this picture — see Ref. ||l^ . An important consequence of 
this model is that the observation of flux jumps depends 
on the stiffness of the vortex structure, i.e. jumps will 
only be observed when collective vortex pinning is strong. 
These flux jumps can thus be used as an indication of the 
existence of a vortex solid of some sort. 

In Fig. 3 we plot close-ups of the (raw) data in Fig. 
2 in the vicinity of 'Birr- Notice that, for small angles 
(Fig. 3a), there is a noticeable kink in the magnetiza- 
tion at an angle dependent fleld, denoted B^ (indicated 
by arrows), just slightly below B^^r- Prior to this kink, 
the flux jumps decay — becoming negligibly small in the 
vicinity of the kink — while, above the kink, the mag- 
netization decreases smoothly to the reversible domain. 
At larger angles (Fig. 3b), flux jumps persist up to the 
kink, but never beyond. In fact, at the largest angles, 
the amplitudes of the flux jumps are very large — much 
larger than the kink amplitude — and, thus, the kink is 
not discernible. However, an abrupt cessation of the flux 
jumps is instead observed at a field whose angle depen- 
dence merges smoothly into the angle dependence of the 
kink field B^ (see Fig. 4a). Thus, we assume that the 
dramatic flux jump cessation, and the kink observed at 
smaller angles, are related. The fact that the flux jumps 
cease at B„i suggests that the stiffness of the vortex sys- 
tem changes — possibly due to a melting transition (see 
below). In what follows, we identify B^ by the kink at 
low angles, and the flux jump cessation at high angles. 

In Fig. 4a, we have plotted B,„ and Bj^r vs. 6. The 
solid lines are flts to the data using the following scaling 
law derived from 2D Ginzburg-Landau theory llO] : 



Bc(6l)sin6i 
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where B^ refers to either B„i or B^rr, and the subscripts 
II and ± refer to the limiting values of these flelds with 
6 — 90° and 6* = 0°, respectively; the same angle depen- 
dence has also been noted for Hc2 10. From the fit to 
Bm(0), we obtain values for Bm± = 3.6 T and B„i|| — 32 
T, at T = 25 mK. Thus, B,„(6l) falls below B^^r, and well 
below /ioHc2 — a region in B,T parameter space where 
no such transition has previously been observed (earlier 
low-temperature experiments measured only Birr 101 )• 
Using two separate models |^,||, Sasaki et al. 10, have 
made a rough estimate of the T = K, 6* = 0° quan- 
tum melting field Bgm± for k-(ET)2Cu(NCS)2. For ei- 
ther model, they obtain Bgjn± ^ 4 T, very close to the 
transition field Bm± we observe at 25mK. Thus, given 
the extremely low temperature of these measurements 
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(T'^Tc/lO^), along with the proximity to the expected 
T = K quantum melting field, we propose that the 
observed transition at Bm(0, T < 200mK) does in fact 
correspond to a quantum melting transition between a 
Q2D vortex lattice phase and a quantum liquid phase. 
The existence of a liquid phase dominated by quantum 
fluctuation effects in the region between B^^r and /ioHc2 
has, indeed, been noted by other authors [ pj| , p4[ . 

The temperature dependence of provides further 
evidence that quantum effects become important as T— s-O 
K. After scaling measurements made at = 47° back to 
9 = 0°, we plot Bm and B^^r as a function of tempera- 
ture in Fig. 4b. In this B,T regime, the irreversibility line 
shows a linear temperature dependence, consistent with 
previous measurements |l6|] , while B^ exhibits a definite 
negative curvature [B„i oc (Tc— T)", a < 1]; no classi- 
cal theory of melting can account for this trend. Never- 
theless, Blatter et al. have observed a similar behavior 
in 2D superconducting films, which they attribute to a 
crossover from quantum melting to a thermally assisted 
dislocation mediated form of melting at higher tempera- 
tures 1^ . The crossover temperature predicted by Blatter 
et al. occurs at around T ~ Tc/lOO, which is precisely 
where we observe the negative curvature (or crossover) 
in the temperature dependence of B,„. Clearly, mea- 
surements spanning a wider temperature range will be 
necessary in order to elucidate the true nature of the 
melting at higher temperatures. However, the fact that 
the phase boundary curves directly towards the T = 
K axis, as T— ^0 K, is strongly suggestive of a quantum 
phase transition. 

Measurements on a second larger sample produced 
similar flux jumps and exactly the same abrupt flux 
jump cessation just below Bi^r, as seen in Fig 3b. In 
fact, a data point for this second sample falls precisely 
on the melting curve in Fig 4a. Although the flux jumps 
are sample dependent [psj , their cessation appears to be 
sample independent, indicating that the supposed melt- 
ing transition is intrinsic, i.e. B„i(T,^) is unaffected by 
sample shape, size, or details of the sample's surface. 

The pronounced hysteresis in B^ could be taken as 
an indication of I''* order behavior, as could the kink 
in M observed at B^; note that, for several traces in 
Fig. 3a, the kink has the appearance of a discontinuous 
jump. However, it is difhcult truly judge the order of the 
transition based on the hysteresis, given the many fac- 
tors which contribute to asymmetry in the M vs. /XqH 
loops. Indeed, these factors may lead to physically dif- 
ferent transitions on up and down sweeps. The approach 
from low fields is preceded by frequent catastrophic 
order flux avalanches and temperature spikes. Thus, dis- 
order and/or dislocations could play a role in the up- 
sweep melting, though the extreme quality of these sam- 
ples would seemingly rule out plastic or glassy behavior 
(our JPR studies have shown the pinning in this material 
to be several orders of magnitude weaker than in the HTS 



JlSf). The approach from high fields, on the other hand, 
simply involves a transition from a weakly pinned liquid 
state to the ordered state. These matters aside, we 
assert that it is the abrupt cessation of flux jumps at B^, 
together with the angle dependence of B^, that suggest 
a 2D melting transition. From the temperature depen- 
dence of B„i, we additionally propose that the melting 
transition becomes I'** order, and is driven by quantum 
fluctuations as T^ 0. 

Taking this new high field/low temperature data, and 
combining it with previous JPR (l3| , local magnetization 
§,|6|, and ^SR |l| measurements, we can construct a 
mixed state phase diagram for k— (ET)2Cu(NCS)2, as 
shown in Fig. 1. It is tempting to connect the present 
data to the data from Refs. |^ and ||l^ with a smooth 
curve (see dashed line in Fig. 1). This would seem to 
suggest that both transitions correspond to Q2D melt- 
ing. We should note, however, that the high tempera- 
ture phase line may be due to a thermally assisted de- 
pinning transition, as noted in Ref. Jl^ . Finally, the 
present measurements have been unable to detect any 
indication of the transition to a Fulde-Ferrell-Larkin- 
Ovchinnikov state reported recently by Singleton et al. 
close to 6* = 90° |6|. 

In conclusion, we have observed magneto-thermal in- 
stabilities in the mixed state of the organic superconduc- 
tor K— (ET)2Cu(NCS)2, which are associated with transi- 
tions between metastable Q2D vortex lattice phases. The 
abrupt cessation of magnetization jumps associated with 
these instabilities serve as an indication of a melting of 
this Q2D vortex lattice phase. Furthermore, this study 
— which is the flrst of its kind in the high-B /low-T limit 
for such a highly anisotropic superconductor — suggests 
that the melting may be driven by quantum rather than 
thermal fluctuations. Future investigations will focus on 
the temperature dependence of this transition, and on its 
evolution away from the high-B/low-T limit. 
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Figure captions 



Fig. 1. The mixed state B,T phase diagram for k— (ET)2Cu(NCS)2. The legends correspond to: <) — ^oHc2 [9]; 
V — 'Birr [16]; • — 2D melting or depinning [13]; □ — first-order transition [3]; A — 3D to 2D crossover [12]; o — 
2D melting (this study). Note that the Q2D sohd (shaded region) and liquid phases occupy most of the available 
B,T— space. 

Fig. 2. Magnetization as a function of magnetic field (up- and down-sweeps) for angles 9 between ~ 30° and 79°; 
the temperature is 25 mK. 

Fig. 3. Magnetization vs. magnetic field at (a) low and (b) high angles 9 (indicated in the figure). In (a), a distinct 
kink may be observed in the magnetization (indicated by arrows) at a field just below Birr', the fiux jumps decay 
smoothly to zero just before B^. At higher angles (b), the kink at B^ is obscured by large amplitude flux jumps. 
However, the abrupt cessation of these flux jumps serves as an alternative indication of the phase transition. 

Fig. 4. (a) Angle dependence of the melting and irreversibility flelds, B,„(0) and Birr (9) respectively; the solid 
lines are fits to the scaling law in Eq. 1. (b) Temperature dependence of the melting and irreversibility fields. 
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FIG. 1. Mola et al. 




FIG. 2. Mola et al. 
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FIG. 3. Mola et al. 




FIG. 4. Mola et al. 



